Strain and structure driven complex magnetic ordering of a CoO overlayer on Ir(lOO) 
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We have investigated the magnetic ordering in the ultrathin c(10x2) CoO (111) film supported on 
Ir(lOO) on the basis of ab-initio calculations. We find a close relationship between the local structural 
properties of the oxide film and the induced magnetic order, leading to alternating ferromagnetically 
and anti-ferromagnetically ordered segments. While the local magnetic order is directly related to 
the geometric position of the Co atoms, the mismatch between the CoO film and the Ir substrate 
leads to a complex long-range order of the oxide. 



PACS numbers: 73.20.-r, 75.70.Ak, 75.70.Rf 

The adsorption of an ultrathin magnetic oxide film, 
such as CoO, on an nonmagnetic substrate offers a rich 
and fascinating playground for studying the interplay of 
the geometrical structure and the magnetic properties. 
Especially if the stable surface orientation of the oxide 
film differs from the substrate, the structural deforma- 
tions of the oxide film and the magnetic order are closely 
related. This relation can be observed in the case of a 
CoO(lll) layer supported on Ir(lOO), where the pseudo- 
hexagonal atomic arrangement of the oxide film has to 
adapt to the square structure of the Ir(lOO) substrate 
surface. Due to the large lattice mismatch of about 9.8% 
between the bulk lattice parameters of CoO(lll) and 
Ir(lOO), epitaxial growth of the CoO overlayer is very 
unfavorable. For the related formation of Coa^O^y [1^ and 
Mn^^O^ [2, 3 films on Pd(lOO) surfaces, this lattice mis- 
match is compensated by the creation of Co vacancies, 
reducing the stoichiometry of the CoO film. Neverthe- 
less, for the adsorption of CoO(lll) on Ir(lOO), the CoO 
film yields the bulk stoichiometry, but is expanded to 
a c(10x2) overlayer, where 9 formula units of CoO are 
supported on 10 square unit cells of Ir(lOO) (Fig[T]). The 
strain-induced relaxations in the CoO overlayer lead to 
sizable structural inhomogeneities, i.e., distorted bond 
angles, different bond lengths, buckling, and different lo- 
cal environments with respect to the substrate geometry. 
Evidently, these structural properties influence both the 
local and long-range magnetic ordering of the CoO over- 
layer, resulting in complex magnetic patterns. 

In the recent years, density functional theory (DFT) 
approaches have evolved as powerful tools to complement 
the experimental efforts to resolve the atomic structure of 
supported surface oxides [4-6 , and to analyze the under- 
lying electronic structure and the concomitant magnetic 
ordering. In this letter, we present extensive DFT calcu- 
lations for two different CoO c(10x2) overlayer structures 
supported on Ir(lOO). In addition, we demonstrate that 
the close relationship between the magnetic configura- 
tions and the structural deformations is already present 
in simplified models of the system, thus illustrating the 



driving mechanism for the magnetic ordering. 

Density functional theory (DFT) calculations were 
performed using the Vienna Ab-initio Simulation Pack- 
age (VASP) [7[ 12 in the projector augmented wave 
(PAW) framework [9]. The exchange correlation func- 
tional was described by the general gradient approxi- 
mation of Perdew-Burke-Ernzerhof (PBE) [10 . To ac- 
count for the localized d-states of the Cobalt atoms, the 
PBE+U approach of Dudarev [11] was used. Varying U- 
J in the range between to 4 eV, we find that a value 
of U-J=l eV leads to the best agreement with the calcu- 
lated geometric structure of the CoO/Ir(100) system with 
experiment. [12 The same value of U-J=l eV was found 
to yield a bulk volume very close to experiment and also 
describes the electronic structure satisfactorily compared 
to a more sophisticated, but very costly, approach using 
hybrid functionals [13 . In addition, we find that the 
most stable magnetic configuration does not depend on 
the choice of U-J. If not stated otherwise all results of the 
present work refer therefore to this choice of U-J=leV. 

A repeated slab model was used consisting of five Ir 
layers and CoO overlayers on both sides of the Ir slab 
and a separating vacuum of 17 A, thus avoiding dipolar 
interaction between the repeated slabs. A Monkhorst- 
Pack type k-mesh of 11x11x1 was used for small cells, 
and a 4x4x1 k-mesh for the c(10x2) structures. In all 
calculations the structural degrees of freedom are fully 
relaxed. 

Recently, an atomistic model has been proposed for the 
c(10x2) CoO(lll) layer on Ir(lOO) surface on the basis of 
LEED experiments [12 . In this study, 59 structural pa- 
rameters had to be optimized to obtain a satisfactory fit, 
demonstrating the complexity of the system, and call- 
ing for corroboration by an up-to-date DFT approach 
such as applied in this work. Most of the experimentally 
derived structural parameters [12] agree with our DFT 
calculations within the experimental uncertainties, with 
only the average deviation of the y [001] coordinates of 
the CoO layer lying outside by 0.01 (O) and 0.02 A(Co) 
^14j . Fig. [ija) presents the result of the fully relaxed 
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FIG. 1. (Color online) Illustration of the calculated structure 
for a c(10x2) CoO(lll) monolayer on Ir(lOO): red: O; blue: 
Co; white: Ir. (a) Experimentally proposed structure; (b) 
competing structure generated by shifting the CoO layer by 
half the Ir-Ir distance along the [010] direction. All directions 
[z,x,y] are given with respect to a (1x1) Iridium cell. See also 
the positional data shown in Fig. [3] 



calculations. In agreement with the experimental data, 
the oxide film is highly corrugated and displays a pro- 
nounced corrugation in terms of height differences of the 
O (-1 A) and of the Co (-0.5 A) atoms. 

In this paper, we also present a second stable over- 
layer geometry, constructed by shifting the overlayer by 
half the substrate Ir-Ir spacing in the [010] x-direction 
(Fig. [TJb)). In this structure, which is of different sym- 
metry, the changes in the height of the O and Co atoms of 
the oxide film are less abrupt, although the corrugations 
are similar to the unshifted case. Surprisingly, the en- 
ergy difference between the these two c(10x2) structures 
favoring the shifted structure (Fig. [TJb)) by 6 meV/Co 
atom is very small and negligible with respect to the 
expected accuracy. Note that a standard PBE calcu- 
lation tilts the balance in favor of the unshifted struc- 
ture (Fig.[lja)) by 8 meV/Co atom. Therefore the DFT 
calculations predict that both phases coexist under suit- 
able experimental conditions. Experimentally, the two 
phases should be distinguishable by scanning tunneling 
microscopy (STM) because the simulated images are dis- 
tinctly different: four bright spots (related to the upper- 
most oxygen atoms) are predicted for the originally pro- 
posed structure [12 , whereas for the new, registry-shifted 
structure five (or three) bright spots should appear. In- 
deed, very recent experimental measurements seem to 
confirm the coexistence of both structures [15 . 




FIG. 2. (Color online) SQE model of a hexagonal (1x1) 
CoO(lll) monolayer on Ir(lOO): red: O; blue: Co; white: Ir. 
(a) Side view along [010]; and (b) top view, (c)-(e) Different 
magnetic configurations; blue: Co majority, light blue: Co 
minority. Configurations (c) and (d) are energetically almost 
degenerate (see text), while (e) is disfavored by 25 meV/Co 
atom. 

The basic driving forces for the structural reconstruc- 
tion can be analyzed in terms of a simplified, quasi- 
epitaxial (SQE) model, where the hexagonal (Ill-like) 
CoO overlayer is accommodated on a square (100) sub- 
strate in a p(lx2) unit cell as shown in Fig. [2] In the en- 
ergetically most favorable structure, the hexagonal CoO 
film binds to the substrate via both the Col and Co2 
atoms and the on-top (01) oxygen atoms at the corners 
of the supercell. Consequently, the Col and Co2 rows 
are pushed together by 0.53 A compared to the ideal po- 
sitions dictated by the substrate, thus forming Co-Ol-Co 
bond angles between 90° (Col-Ol-Col) and 135° (Col- 
01-Co2), that resemble the quasi-hexagonal arrangement 
reminiscent of hex-BN networks. The second type of oxy- 
gen atoms (02), located in bridge sites, is tilted out of the 
plane and hence does not directly contribute to the bond- 
ing to the surface. The oxygen atoms display a 01-02 
buckling of 0.97 A, forming a three-sided pyramid with 
Co-O-Co bond angles close to 90° such as found in the 
bulk rock-salt structure along the [111] diagonal. Such a 
buckling is much less pronounced for the Co atoms, which 
are either located in hollow sites (Col), with a vertical 
distance of 1.94 A above the substrate, or in bridge sites 
(Co2), at a height of 2.27 A (cf., Fig.[2[a)). For compari- 
son, a pure metallic Co overlayer on Ir(lOO) is also found 
to occupy hollow substrate sites at a height of 1.68 A. 

The different coordination of the Co atoms is directly 



3 



reflected in the magnetic ordering: The spin moments 
of the Col atoms in the hohow sites ahgn ferromagneti- 
cahy, comparable to the ferromagnetic order of the pure 
metahic Co overlayer (Fig. [2jc), dark blue atoms). On 
the other hand, the magnetic order of the weaker bonded 
Co2 rows on the substrate bridge sites is less well de- 
flned: the difference in energy between the antiferromag- 
netic ordering and the ferromagnetic ordering within the 
Co2 rows — accompanied by an antiferromagnetic cou- 
pling to the Col rows — is only 0.2 meV per Co atom 
(Figs.[2l;c,d)). 

All other reasonable (collinear) magnetic orderings are 
unfavorable by at least 25 meV per Co atom, including 
antiferromagnetic alignments within the Col rows or fer- 
romagnetic coupling of all Co atoms. The induced polar- 
ization on the nearby Ir atoms is small since the neighbor- 
ing Co atoms have partially opposite magnetic moments. 
The same holds for the magnetism induced in the oxy- 
gen atoms. Therefore we conclude that the Co atoms 
located on hollow sites have a strong tendency for fer- 
romagnetic ordering along the chains, and antiferromag- 
netic coupling across the lines of O atoms to the bridge- 
site Co2 chain. Applying the Kanamori-Goodenough- 
Anderson p^BHTS] rules in a simple way, one would pre- 
dict a ferromagnetic coupling of the Col and Co2 chains 
connected by 02 atoms via direct Col-Co2 overlap (Col- 
02-Co2 bond angles close to 90° and Col-Co2 distance 
of 2.62 A), while the coupling involving the planar 01 
type atoms is less clear: The 90° Co2-01-Co2 angles and 
a Co2-Co2 distance of 2.74 A would again favor direct 
overlap ferromagnetic ordering along the Co2 rows, but 
the 135° Col-01-Co2 angles and Col-Co2 distances of 
3.56 A would certainly work against ferromagnetic cou- 
pling, although antiferromagnetic coupling involving the 
01 atom would be strongest for 180° angles. However 
such a simple picture for the CoO fllm is certainly com- 
plicated by the interaction with the substrate, either by 
direct coupling of the Co spins via the substrate, or by 
changes in the oxygen orbitals caused by bonding to the 
substrate that would modify the super-exchange interac- 
tion. Thus, the magnetic ordering observed for the CoO 
fllms should be characterized as a competition between 
ferromagnetic coupling between the Co atoms by direct 
overlap (and also mediated by the substrate) and antifer- 
romagnetic coupling across oxygen atoms in the hex-BN 
like structures of the layer. 

However, the small energy differences are of the same 
order as contributions from non-collinear spin arrange- 
ments. Since the full c(10x2) structure is too complex 
for a non-collinear study we assess this possibility by per- 
forming calculations for selected non-collinear orderings 
between Col and Co2 spins based on the SQE model 
(Fig. 2). Considering Fig. 2c, as a flrst step (corre- 
sponding to the collinear spin-polarized calculation) the 
collinear spin directions are flxed along [100] (out-of- 
plane) and then in a second step the FM Col or the 



c(10x2)CoO/lr(100) 




^ 1.8 
°< 1.6 
>?1.4 

3.1 

n"2.3 

1.9 
_ 0.3 

•i„o.o 

o 

^-0.3 
^ 1.8 
"S1.6 
1.4 



02(SQE) 




Co 1 (SQE) 



Co1(SbE) 




-01,2(SQE) --- 



(C) 



(d) 



(e) 



(f) 



-10 -5 5 

X[010] (A) 



10 



FIG. 3. Structural details of a hexagonal (1x1) CoO (111) 
monolayer on Ir( 100). (a-c): Comparison between experiment 
and PBE+U (U-J=leV) results, (d-f): CoO layer shifted by 
half of the substrate Ir-Ir spacing along x [010]. Thin vertical 
lines show the ideal lateral positions of the substrate Ir atoms 
along X [010], and empty (filled) diamonds the ideal hollow 
(top) sites of the Ir substrate. For comparison, the "SQE" 
labelled dashed horizontal lines denote the atomic positions 
for the SQE model of Fig. [2] The z [100] (height) values 
for CoO in (a) and (d) are given with respect to the mean 
Ir positions in the interface layer, while in (b-c) and (e-f) 
the y [001] values are given with respect to the ideal hollow 
sites of the Ir substrate where Co atoms of a metallic Co 
overlayer would be found. To allow for a better comparison 
with experiment the calculated lateral distances have been 
scaled to the experimental Ir lattice constant of 3.84 A (theory: 
3.88 A). 



AF Co2 spins are rotated by 90° perpendicular to [010] 
(in-plane). The energy gain for these non-collinear con- 
flgurations is about 8 meV/Co. Estimating the coupling 
to the substrate one flnds that exchanging the spin axis 
of Col and Co2 favors a Col in-plane configuration, but 
only marginally by 0.6meV/Co. Rotating the Co2 spins 
in-plane, while keeping the FM order (Fig. 2d), desta- 
bilizes the configuration by 20 meV/Co with respect to 
the one of Fig. 2c. Therefore it seems likely that non- 
collinear spin ordering occurs in the full c(10x2) struc- 
ture, but that the magnetic coupling to the substrate 
plays only a minor role here. 

The main patterns seen in the SQE model also hold for 
both of the much more complex cases of c(10x2) CoO 
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FIG. 4. (Color online) Magnetic structure of a c(10x2) 
CoO(lll) monolayer on Ir(lOO): red: O; blue: Co majority; 
light blue: Co minority; white: Ir. The top panel shows 
the experimentally proposed c(10x2) structure [121 114] , the 
bottom panel a competing structure with a registry-shifted 
c(10x2) CoO overlayer, which is slightly favored by 6 meV/Co 
atom. 

overlayer structures on Ir(lOO) (see Fig. [4|. Although 
the mismatch between the CoO film and the Ir(lOO) sub- 
strate leads to variations in the local configuration, the 
segments of the Co rows, either positioned on bridge-like 
or on hollow-like sites, are clearly visible in the fully re- 
laxed structure (Figj3|. However, as the c(10x2) CoO 
film is strained in the [010] x-direction compared to the 
SQE model, there is a continuous transition from the 
hollow-like to the bridge-like Co segments, and hence the 
Co atoms are not located in the ideal positions of the 
SQE model. Consequently, the chains in the [010] x- 
directions get broken, when Col atoms are shifted from 
their ideal hollow positions to bridge positions; similar 
arguments holds for the Co2 chains. Therefore, the full 
model does not display infinite chains of one type of Co 
atoms along the [010] x-direction, but rather alternat- 
ing segments of both types arrange in a zig-zag pattern 
(compare Figs, [s] and [4|. Still, the height and in-plane 
modulations for both the O and Co atoms shown in Figs, 
[sja-c) match quite nicely the coordinates of the corre- 
sponding atoms of the SQE model which already captures 
the deviations from the ideal positions perpendicular to 
the Co/0 chains ([001] y-direction) . Also the agreement 
with the experimental LEED results [121 [14] is excellent. 

The predicted correlation between the spin alignment 
and the local geometry is also observed for the full model. 
For both CoO overlayer structures, the calculations show 
a particularly pronounced binding to the Ir substrate if 
the oxygen atoms are in on-top sites and the neighbor- 
ing Co atoms are in hollow sites, forming anchor sites 



around which the less binding Co2 chain fractions are 
placed. This leads to the formation of two types of Co 
segments, a lower segment consisting of 4 (Col) atoms 
in hollow positions, and a higher segment of 5 (Co2) 
atoms. As already observed in the SQE model, all Col- 
like atoms are ferromagnetically aligned (Fig. [4j upper 
panel). Compared with the SQE model, we also pre- 
dict antiparallel ordering of Co2 (bridge-like) atoms with 
an enhanced energy difference between anti- and parallel 
ordering within the same Co2 segment of 6 meV per Co 
atom. Therefore it is likely that the atoms connecting the 
two different chain segments (Col-like and Co2-like) have 
a non-collinear alignment. A very similar relation can be 
observed for the second, registry-shifted CoO structure 
(Fig. [4j lower panel). Due to the shift in the substrate 
lattice, there are now segments with an odd number of 
Col atoms (now 5) and an even number of Co2 atoms 
(now 4), and vice versa for the O atoms. For this geom- 
etry, the spin alignments are just the ones expected from 
the SQE model, i.e., all orderings within the chain seg- 
ments are now parallel. As the non-collinear calculations 
for the SQE model show, a rotation of a ferromagnetic 
Co2 spin segment is unfavorable. Hence non-collinear 
effects are expected to play a lesser role for the registry- 
shifted phase. Concerning the chain segments of the 5 
high lying O atoms, the border atoms of these chains are 
lower in height because they are significantly shifted lat- 
erally towards the on top substrate sites. Overall, this 
leads to a surface of smoother corrugation (cf., Figs.[TJb) 
and [sj^d-f)) compared to the unshifted case. Yet, the 
small energy difference between the plain registry-shifted 
c(10x2) phases can lead to the coexistence of both struc- 
tures, and hence incommensurate spin ordering. 

In summary, we have investigated and analyzed the 
complex magnetic ordering in the supported c(10x2) 
CoO films on Ir(lOO). The calculations predict two 
closely related CoO overlayer phases, built up by Co seg- 
ments consisting of either 4 or 5 atoms. We find a close 
relationship between the structural relaxations and the 
local magnetic ordering in the overlayer: the segments of 
Co atoms above hollow-like positions of Ir(lOO) clearly 
favor ferromagnetic coupling along the lines and anti- 
ferromagnetic coupling between the ferromagnetic rows, 
while the magnetic order is less pronounced in the seg- 
ments with bridge-like atoms. Based on model calcula- 
tions, non-collinear arrangements are expected to occur 
predominantly for the AF ordered chain segments. On 
the basis of these findings and their analysis, we expect a 
similar relationship between structure and magnetic or- 
dering for related transition-metal oxide overlayers. Con- 
sequently, inducing a specific deformation by the cover- 
age by the oxide overlayer may provide a strategy for 
designing magnetic properties of complex surface oxides. 
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